Paramagnetic reentrance of the ac-screening: Evidence of vortex avalanches in Pb 

thin films 
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We have studied the influence of a square array of pinning centers on the dynamics of vortex 
avalanches in Pb thin films by means of ac- and dc- magnetization measurements. Close to the su- 
perconducting transition T c the commensurability between the vortex lattice and the pinning array 
leads to the well known local increments of the critical current. As temperature T decreases, match- 
ing features progressively fade out and eventually disappear. Further down in temperature vortex 
avalanches develop and dominate the magnetic response. These avalanches manifest themselves 
as jumps in the dc-magnetization and produce a lower ac-shielding giving rise to a paramagnetic 
reentrance in the ac-screening \'{T). Within the flux jump regime two subregimes can be identified. 
Close to the boundary where vortex avalanches develope, the field separation between consecutive 
jumps follow the periodicity of the pinning array and a field and temperature dependent screening 
is observed. In this regime, the response also depends on frequency / in agreement with theoretical 
models for magnetothermal instabilities. At low enough temperatures and fields, the screening sat- 
urates to a constant value independent of T, H , and /, where jumps are randomly distributed. We 
have also found that vortex instabilities occupy a larger portion of the H — T diagram in patterned 
samples than in films without nanoengineered pinning sites. Finally, we discuss the possible origin 
of the vortex avalanches and compare our results with previous experimental and theoretical studies. 



C 

O 
o 



> 

CO 

o 
o 

-I— > 

C3 



I 

G 

O 

o 



X 



I. INTRODUCTION 

If a type-II superconductor is cooled down in a zero 
applied field (ZFC experiment) and subsequently an ex- 
ternal field H larger than the first critical field H c \ is 
applied, flux-bearing vortices enter through the sample's 
borders until they are captured by pinning centers. As a 
consequence, the system achieves an inhomogeneous flux 
distribution with a higher density of vortices near the 
borders that progressively decreases toward the center of 
the sample. The spatial variation of the locally averaged 
field B(r) gives rise to supercurrents J in the sample that, 
in the stationary state, accommodate to be exactly the 
critical current J c everywhere. The resultant inhomo- 
geneous flux distribution of this so-called critical state, 
represents a self organized state that under small pertur- 
bations (like local temperature fluctuations <5Ti) can lead 
to vortex avalanches in order to maintain J = J C {H,T) 
in that region. Because of the dissipation produced when 
flux lines move, each avalanche gives rise to a heat pulse 
and therefore a local increment of the temperature STf. 
If 6Tf < STi the critical state remains stable under these 
perturbations, otherwise avalanches of vortices draw the 
sample to a highly resistive state, 1 

The way in which the vortex lattice reacts to the local 
overheating is crucial to determine the subsequent dy- 
namics of the system. If the diffusivity of the magnetic 
flux (given by the resistivity p) is smaller than the ther- 
mal diffusion coefficient (given by k/C, where K the heat 
conductivity and C is the heat capacity), the hot spot 
propagates in a frozen magnetic and current distribution 
(dynamic approximation). On the other hand, if ther- 
mal diffusivity is smaller than magnetic diffusivity there 



is not enough time to remove or distribute the heat pro- 
duced by the vortex motion (adiabatic heating) ? In both 
cases, the stability criterion for the critical state reads, 
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where s is the characteristic sample dimension, e = C/fM) 
in the adiabatic approximation and e = n/p in the dy- 
namic approximation^ 

At high temperatures and fields, where the critical cur- 
rent is small, vortex avalanches are usually not seen. 
As temperature and field decrease, critical current in- 
creases and below a certain boundary H*(T) vortex 
avalanches develop. In principle this scenario should ap- 
ply with minor differences for both bulk samples and thin 
films. However, recent magneto-optical images (MO) 
have shown that in thin films, with H applied perpen- 
dicular to the plane of the film, the flux pattern in the 
sample exhibits a richer morphology than the smooth 
progressive flux penetration observed in the bulk. For 
this particular geometry, flux invasion occurs via den- 
dritic structures which cover a substantial portion of the 
sample's area and grow very fast (i.e. under adiabatic 
conditions)^ The thermal origin of the observed insta- 
bilities has been recently demonstrated by strongly sup- 
pressing the dendrite instabilities as the thermal contact 
is improved^ 

Typically, the observed sudden penetrations of the flux 
front into the sample are accompanied by sharp jumps in 
the dc-magnetizationi£<£<£ These jumps generate a noisy 
response which can undermine the technological applica- 
bility and perspectives of superconducting devices at low 
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temperatures. However, recent promising experimental 
results have shown that this noise can be substantially 
reduced by introducing an array of pinning centers jSaLifl 
at expenses of increasing considerably the region in the 
H — T diagram where flux jumps occur 

Additionally, in a recent work, Aranson et al. 12 have 
theoretically predicted that the presence of a periodic 
modulation of the critical current would give rise to a 
growth of the branching process of the dendrites in com- 
parison with a homogeneous pinning distribution. Ev- 
idence supporting this picture was reported by Vlasko- 
Vlasov et aL 13 who analyzed the MO images of Nb films 
patterned with a square lattice of holes. The authors 
show that as the field is progressively ramped up, first 
the flux enters from the edges in stripes with boundaries 
along the principal axes of the pinning array and then 
new stripes jump between previously developed stripes. 
The width of these stripes involves several units cell of 
the pinning array. 

On top of that, magnetization measurements per- 
formed on Nbi^ and Pbil films with a square array of 
holes show that the field separation between consecutive 
jumps commensurate with the period of the underlying 
pinning array. This is a surprising result since matching 
features are typically seen only very close to T c where 
intrinsic pinning and self- field effects are not relevant. 
Clearly, the influence of the vortex pinning in the mor- 
phology and dynamics of the flux penetration on these 
kind of systems is an issue that has not been fully ad- 
dressed so far and deserves further investigations. 

In this work we study the flux-jump regime in Pb thin 
films with and without periodic pinning by means of ac- 
and dc- susceptibility measurements. In the samples with 
a square antidot array several regimes can be identified 
as a function of temperature. For temperatures T < T c , 
matching features appear when the vortex lattice com- 
mensurates with the pinning array. As temperature de- 
creases these effects progressively fade out. At a certain 
field-dependent temperature the shielding power of the 
sample is dramatically reduced and a reentrance in the 
ac-screening x'(^) is observed. This effect is accompa- 
nied by a substantial increase of the dissipation \" ■ We 
demonstrate that the observed reentrance is related to 
the appearance of flux jumps in the sample. Finally we 
discuss the origin of the quasi-periodic jumps in terms of 
a multiterrace critical state. 



Si02 substrate which allows us to perform a direct and 
reliable comparison in order to discriminate the effects 
of the pinning array. From the temperature dependence 
of the upper critical field H C 2(T) we have estimated a 
superconducting coherence length £(0) = 33 ± 3 nm for 
all samples. 



TABLE I: Lateral dimensions (wi and W2), thickness (t), and 
critical temperatures (T c ) for all the films studied. AD stems 
from square array of antidots and BH from square array of 
blind holes. 



Sample 


wi (mm) 


W2 (mm) 


t (nm) 


Tc (K) 


AD 15 


1.9 


2.0 


13.5 


7.10 


AD65 


2.3 


2.5 


65 


7.21 


BH100 


3.2 


3.3 


100 


7.22 



The details of the sample preparation can be 
found in Ref.0]. Briefly, the predefined resist- 
dot patterns were prepared by electron-beam lithogra- 
phy in a polymcthyl metacry late/methyl metacrylate 
(PMMA/MMA) resist bilayer covering the Si0 2 sub- 
strate. A Ge(20 A)/Pb/Ge(200 A) film was then 
electron-beam evaporated onto this mask while keeping 
the substrate at liquid nitrogen temperature. Finally, 
the resist was removed in a lift-off procedure in warm 
acetone. The BH100 blind hole array was fabricated by 
depositing an additional 25 nm-thick Pb film on top of a 
75 nm-thick Pb film with an array of antidots. 

The ac-susceptibility measurements x(H, T) — x'+*x" 
were carried out in a commercial Quantum Design-PPMS 
device with drive field amplitudes h ranging from 2 mOe 
to 10 Oe, and the frequency / from 10 Hz to 15 kHz. 
The data were normalized to have a total step A\' = 
1, with H = at low temperatures and ac drives. This 
system provides a temperature stability better than 0.5 
mK which is crucial for measurements near the criti- 
cal temperature. The dc-magnetization measurements 
were obtained using a QD-MPMS SQUID magnetometer 
equipped with a 5 T magnet. 



III. RESULTS AND DISCUSSION 



A. Dc-magnetization 



II. EXPERIMENTAL DETAILS 

The experiments were conducted on Pb thin films with 
different pinning arrays. The dimensions and critical 
temperature for each sample are summarized in Table 
In all patterned samples the square antidot array 
consists of square pinning sites with lateral dimension 
b = 0.8 /im and period d = 1.5 /im which corresponds to 
a first matching field H\ = 9.2 G. Simultaneously with 
each patterned film we deposited a plain reference film on 



In order to identify the temperature range where flux 
jumps occur, we have measured zero field cooled (ZFC) 
and field-cooled (FC) dc-magnetization for a fixed field 
H. The result of these measurements is shown in the 
main panel of Fig. \T\for the AD15 sample at H = 5 G. 
A reversible response is obtained for T c — T < 0.7 K, a 
temperature range higher than the expected for the irre- 
versible line at this field^ probably due to the presence 
of an undesirable remanent field in the initial cooling pro- 
cedure. 
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FIG. 1: Main panel: Zero-field cooling (filled circles) and 
field-cooling (open circles) dc-magnetization as a function of 
temperature for the AD 15 sample at H = 5 G. The vertical 
black arrows indicate the temperature where the loops shown 
in the insets where recorded. 



The AD 15 sample is the same as that used in our pre- 
vious studies^ In that work we show that commensura- 
bility effects between the flux line lattice and the pinning 
array manifest themselves as a peak in the critical current 
or as a higher screening in the ac-response. In this partic- 
ular sample, well defined matching features appear only 
for T c -T < 0.6 K. For lower temperatures (T < 6.5 K), 
the decrease of the penetration depth A(T), the growth of 
the intrinsic pinning strength and self-field effects^ lead 
to a less pronounced matching features, which eventually 
disappear. Decreasing further the temperature there is a 
clear transition to a regime where magnetization becomes 
notably noisy This transition occurs for both ZFC and 
FC curves around T — 5.7 K indicating that the onset 
of the crossover to a smooth behavior at high tempera- 
tures is highly reproducible. Of course the noisy response 
(8m ~ 0.1m) is not related to the experimental resolution 
of the used device since in this regime the signal/noise 
ratio is higher than near T c where a smooth curve is ob- 
tained. Instead, we argue that at T — 5.7 K a transition 
to a more dissipative state owing to flux jumps takes 
place. Evidence corroborating this interpretation is ob- 
tained by measuring hysteresis loops at T = 4.8 K deep 
in the noisy regime (upper inset of Fig. ^| and at T = 6 
K where the noise is absent (lower inset of Fig. It can 
be seen that at T = 6 K jumps in the dc-magnetization 
never occur. In contrast to that, at T = 4.8 K flux 
jumps are present in the low field region up to ~ 100 G, 
where the average magnetization reaches a maximum. 
For H > 100 G, the irreversible magnetization decreases 
smoothly as field increases in agreement with a critical 
state scenario. As it has been pointed out previously, 



the magnetization peak at 100 G indicates the onset of 
vortex avalanches in the sample. 6 



B. Ac-susceptibility 

The appearance of flux-jumps should also be reflected 
as a lower efficiency to screen out an external ac field. 
This effect can be seen in the lower panel of Figure El 
where the screening x' and the dissipation x" a s a func- 
tion of temperature are shown for the AD15 sample at 
H = 5 G and several ac excitations. 




FIG. 2: Temperature dependence of the ac-susceptibility \ = 
x' + ix" a t several ac-drives h for a plain film (upper panel) 
and a film with a square array of holes (lower panel). In 
the lower panel the transition to a flux-jumps regime appears 
as a reentrance signaled by the features Tdvp and T pea k (see 
arrows) . 

At low ac drives no features indicating the transition 
to a flux-jumps regime are observed. For /iMGa para- 
magnetic reentrance in the screening at T ~ 6.3 K signals 
the onset of flux-jumps regime. It should be noted that 
the ac drive at which the reentrance in x'(T) first ap- 
pears, h ~ 1 G, is the same order as the average distance 
between the jumps Hj ~ 3 G. This result suggests that 
as long as h <C Hj new vortex avalanches are not trig- 
gered and the transition to the flux jumps regime is not 
detectable. Consistently, we have observed a weak influ- 
ence of the amplitude h on the position where jumps first 
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develop. This is shown in the inset of Figure |21 where we 
can see that at low amplitudes, both the local minimum 
Tdip and the local maximum T pea k, remain almost con- 
stant as h increases. From h = 4 G up, the transition 
temperature slowly decreases with increasing h. As we 
will show below this effect is a consequence of a Tdip de- 
creasing with increasing the total field H + h. There is 
a systematic small discrepancy between the onset of the 
flux-jumps regime determined by dc-magnetization and 
ac-susceptibility. This effect can be attributed to a faster 
onset of vortex avalanches due to the ac shaking. 

It is worth to note that the temperature evolution of 
the curves at temperatures above Tdip and below T pea k 
follows a quite different trend. Similar behavior has been 
recently reported by Passos et ali& for MgB2 samples and 
attributed to a possible field induced granularity. The 
origin of the reentrant behavior in the screening can be 
understood from the upper inset of Fig. ^ In this figure 
we can see that the development of avalanches leads to 
a reduction of the average critical current and therefore 
also to a reduction in the screening properties. In other 
words, when the system crosses from a "non-jumpy" 
regime to a "jumpy" regime, either by sweeping field or 
temperature, a strong suppression of the screening power 
occurs. As a rule, every peak in the dc-magnetization (ei- 
ther in field or temperature) will manifest itself as a dip 
in the ac-screeningi£2i2i 

For comparison, in the upper panel of Figure [21 we 
show similar measurements performed on a plain film 
(without nanoengineered pinning array) for the same ac 
drives and dc field. In the window of temperatures shown 
here, no features indicating the flux-jump transition are 
found. However, a reentrance appears at Tdip ~ 3.8 K for 
H = 5 G. This reduction of Tdip in the unpatterned sam- 
ple is in agreement with recent magnetization measure- 
ments performed in samples with and without a periodic 
pinning array. 11 

It should be noted that performing ac-susceptibility 
measurements the "noise" produced by the flux jumps 
is removed. This is so because the measurements here 
presented were obtained at a relatively high frequency 
/ = 3837 Hz and with an integration time of 1 sec, there- 
fore the resultant \ reflect an average after cycling 3837 
times a minor loop. By reducing the number of cycles 
one may approach to the case of dc-magnetization mea- 
surements and jumps become more apparent. 

The transition to the flux-jump regime becomes more 
obvious by plotting the imaginary versus the real part of 
the ac-susceptibility: \" x x' as shown in the main panel 
of Fig. [21 The advantages of this so-called Cole-Cole plot 
is two fold, first it makes possible a comparison among 
different samples without knowing the specific value of 
the critical current^ secondly it allows one to identify 
clearly different vortex dynamics regimes. In Figure |3| 
we show this representation with x obtained by sweeping 
temperature at several h (part of these data were already 
shown in Figure |5J). All these curves show a similar gen- 
eral behavior. Starting from high temperatures (right 




FIG. 3: Main panel: \" vs \' obtained from x(T) curves at 
several ac-amplitudes h for the AD 15 sample (some of these 
curves are shown in the lower panel of Figure |2j. The small 
arrows indicate the evolution of the curve as T decreases. The 
inset shows the transition temperatures Tdip and T pea ^ (see 
Figure |5J as a function of h. 



side) we first observe an increase in both dissipation and 
screening as T decreases, until the transition temperature 
Tdip is reached, then a sudden reduction of the screening 
together with an increase of the dissipation occurs down 
to T pea k. Below this temperature, the curves continue 
with a smooth evolution. Remarkably, the whole family 
of curves merges in two well distinguished evolvents, one 
for T > Tdip and the other for T < T pea k , the latter being 
more dissipative than the former. Additionally, the peak 
of maximum dissipation is shifted from \' ~ —0.365 for 
T > T dlp to x' ~ -0.408 for T < T peak . This increase 
in magnitude of the maximum x" together with the shift 
towards x' — — 1 is commonly observed in high tempera- 
ture superconductors and attributed to creep effects. As 
we will discuss below the higher dissipative state is ac- 
companied by a /-dependence susceptibility analogous 
to a creep regime where the frequency dependence of the 
magnetic response appears as a consequence of the dy- 
namic evolution of the system (not considered in a con- 
ventional critical state model). 

Once we have identified the temperature Tdip (or 
Tpeak) as the transition temperature to a flux-jumps 
regime, we can determine the portion of the H — T dia- 
gram where these instabilities dominate by tracking the 
dip position as a function of the dc field. Some of these 
curves are shown in the main panel of Figure 0] for the 
AD65 sample at h = 3 G. In this figure it can be seen 
that the transition Tdip shifts to lower temperatures as 
H increases. This result is a consequence of a critical 
current J c decreasing with increasing H and in agree- 
ment with previous experimental reports^i It is also in- 
teresting to notice that at low temperatures and fields 
all these curves merge into a single universal behavior. 
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This field independent curve indicates that the average 
irreversible magnetization in this region is also field inde- 
pendent. This is confirmed by the dc magnetization mea- 
surements shown in the main panel of Figure [S] for the 
BH100 sample and by recent experiments performed on 
plain Pb filmsi^i Another important observation is that 
at low temperatures the universal magnetization seems 
to saturate to a value x' ~ —0.91, and hence it also 
becomes temperature independent. Of course, at higher 
fields (H > 100 G) this field independent behavior is pro- 
gressively lost as the system approaches the boundary 
H*(T). 

This striking T- and i/-indepcndent average magneti- 
zation at low temperatures might be attributed to the 
coexistence of two different species of vortices. Indeed, 
it has been shown that in this regime dendrites have a 
finger-like structure with a size which remains unaltered 
with further changing the field. This very stable pattern 
of flux penetrated regions is a consequence of the ten- 
dency of dendrites to avoid each other. Once the field 
is reduced, already developed dendrites may help to re- 
move flux from the sample and also channel incoming 
vortices when the field is again increased. In this pic- 
ture, the ac response would be dominated by the easy 
motion of vortices inside the dendrites thus leading to a 
lower saturation value \' ~ —0.91. 
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FIG. 4: Main panel: x'{ T ) for tne AD65 sample at h = 3 
G and several fields. The vertical dashed line indicates the 
lowest temperature where matching features are still barely 
defined. The inset shows the different observed regimes in 
a H — T phase diagram for the AD15 (square symbols) and 
AD65 (circle symbols) samples. The boundary H*{T) repre- 
sents the position of the local minimum in the x'(T) curves 
shown in the main panel. 

In the inset of Fig. 0] we show the resultant phase di- 
agram obtained following the procedure described above 
for the AD65 (open circles) and AD15 (filled squares) 
samples together with the upper critical field H C 2 (T) for 



both samples^ An almost linear temperature depen- 
dence of the boundary Hd ip (T) is observed in the whole 
range of temperatures. Note that the observed different 
slopes of the boundary H * (T) (a factor of 3 steeper for 
the AD65 sample) cannot be ascribed to a different tem- 
perature dependence of H C 2 since the upper critical field 
is roughly sample independent. As a consequence, the 
H — T region of flux-instabilities for the AD 15 turns out 
to be smaller than for the AD65 sample. This reduction 
of the flux-instability region with decreasing the film's 
thickness is in agreement with previous observations in 
Nb thin films. 6 Indeed, if H is applied perpendicular to 
the sample's surface, the effective size s which determines 
the stability criterion in eq.lfj^) can be approximated as 
s ~ \J wt/2, where w is the width of the film&2^ Us- 
ing the sample dimensions shown in Table [3 we obtain 
s = 3.8 pm and s — 8.9 ^m for the AD15 and AD65 
samples, respectively. From eq.lfjQ) we have that for di- 
mensions s such that s 2 > s 2 crit — IfC-gjH vortex instabil- 
ities appear. In other words, the smaller the s the lower 
the boundary H* (T) , in agreement with our observation. 

As we pointed out previously, the transition to the 
vortex-instability regime manifests itself as a peak in the 
dc magnetization loops such as that shown in the upper 
inset of Fig.^ Similarly, we can determine this transition 
field by performing x'(-ff) measurements at several tem- 
peratures as is shown in the inset of Fig. for the AD65 
sample. In this figure, the data recorded at T = 6.5 K 
fall, for all fields, in the stable regime where no jumps 
are observed and therefore exhibits the standard maxi- 
mum screening at zero field. For the rest of the explored 
temperatures T < 6 K the maximum screening is no 
longer located at zero field but at H pea k which roughly 
coincides with the boundary H*(T) shown in the inset of 
Fig. We observe that at low temperatures a plateau 
in x'(H) appears around H = in agreement with our 
previous remark. This constancy of x' {H) implies a fixed 
size of the flux avalanches and a field and temperature 
independent magnetization. This is consistent with re- 
cently reported avalanche distribution performed in sim- 
ilar samples^ and magnetization measurements on plain 
Pb films^ According to Ref-O, for T < 5 K the size 
distribution of flux-jumps does not depend strongly on 
temperature and in particular, the maximum jump re- 
mains almost constant. This effect is not inherent of 
samples with a periodic pinning array since similar re- 
sults were found in unpatterned Pb films^ At lower 
temperatures (T < 3 K) a peak in x'(H) around H = 
reappears although its maximum screening at H = is 
always smaller than x' = —0.91. So, still this saturation 
value imposes an upper bound for the maximum possible 
screening in the flux-jump regime. 

Lets now discuss what kind of instability corresponds 
to the observed behavior. An estimation of the ratio 
t = na/C between thermal and magnetic diffusivities, 
~ C = 0.67 J/Kg-K and k = 3.4 W/cm-K, and 
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considering that p ~ pff where pff = p n H/H C 2 is the 
flux-flow resistivity and p n ~ 0.02 pflm the normal state 
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FIG. 5: Main panel: x'(T) for the AD65 sample at zero field, 
h — 3 G and two extreme frequencies. The inset shows the 
field dependence of the screening \' f° r several temperatures. 
At low T and H, the flux-jump regime manifests itself as a 
reduction of the screening. 



resistivity, gives r ~ 1 at 4 K and H = 100 G. There- 
fore, within this context neither the adiabatic nor the 
dynamic (or isothermal) approximations seems to be ap- 
propriate. However, in the estimate we made several im- 
portant points are not included. First, the heat removal 
through the substrate which results in a lower local in- 
crease of the temperature^ Second, the field and tem- 
perature dependence of the specific heat and the thermal 
conductivity. All these effects tend to reduce the value of 
the effective r thus approaching to the adiabatic limit. In 
addition, our estimate is based on a critical state model 
where the dendritic penetration observed in these sam- 
ples is not considered^ Moreover, it is important to note 
that flux jumps are not mounted on a continuous prolon- 
gation of the curve observed at high field (see Fig@ , thus 
suggesting that the observed vortex avalanches cannot be 
described within a critical state scenario, ft is believed 
that in thin films the long-range vortex- vortex interaction 
and the non-local current-field relation are key ingredi- 
ents which should eventually be incorporated in order to 
describe the branching process of the flux penetration! 3 ^ 

On the other hand, the similarity between the mea- 
surements presented in this work and those previously re- 
ported in several other materials (MgB2, 7 Nb^i 3 . Pb^ 3 . 
Nb 3 Sn22) points to magnetothermal instabilities as the 
origin of the flux jumps at low temperatures. As we men- 
tioned in the Introduction, there is a rich zoology of pos- 
sible vortex avalanches like smooth flux fronts, finger-like 
penetration or highly-branched tree-like flux invasion. In 
MgB2, MO images indicate that quasi- unidimcnsional 
fingers occur at low temperatures whereas at intermedi- 
ate temperatures a branched structure dominates. Simi- 
larly, Vlasko- Vlasov et ah — showed that in Nb films pat- 



terned with an array of holes flux penetrates as stripes 
at low fields (resembling the fingers observed in MgB2) 
whereas at higher fields a branching process produced 
by "magnetic discharging" connects the original stripes. 
This fast flux penetration strongly suggests that our re- 
sults can be described within the adiabatic approxima- 
tion. 

Another interesting result comes from the analysis of 
the frequency dependence of the screening. According 
to the commonly used theoretical models, in both limits, 
dynamic (t>1) and adiabatic (r <C 1), the size and the 
distance between jumps should decrease with increasing 
the field sweep rate dH/dtm^ It has been also experi- 
mentally shownSi that the mean value of the magneti- 
zation, which is in fact what we are able to determine 
with ac-susceptibility, decreases with increasing dH/dt. 
In order to study this effect we have measured the x' (T) 
at h = 3 G for two extreme frequencies / = 56 Hz and 
/ = 10 kHz, as shown in the main panel of Fig. [3] This 
frequencies correspond to a field sweeping rate of 0.05 
T/s an 12 T/s, respectively, thus widely covering the 
range where /-dependence has been observed in other 
materials^ Surprisingly, there is almost no frequency 
dependence unless in a narrow temperature window right 
below Tdip where a smaller screening is detected for the 
higher frequency as expected^ On the other hand, the 
lack of dH/dt dependence in the low temperature regime 
has been also reported in Ref. |MTT| although in these 
publications the used dH/dt covered a much smaller 
range. 



C. Influence of the periodic pinning 

We now turn to the analysis of a possible matching be- 
tween the periodicity of the jumps and that imposed by 
the pinning landscape. It has been recently shown, first 
by Terentiev et al*^ in Nb films and later on by Hebert 
et a/^i in Pb films that within the flux-jump regime and 
at high enough temperatures and fields, the distance be- 
tween consecutive jumps Hj coincides with a multiple of 
the matching field nH±, with n integer. We have con- 
firmed that this effect seems to be also present in the 
BHI00 sample. Indeed, the upper right inset of Fig. 
shows a zoom-in of the decreasing branch in the flux- 
jumps region for the BH100 sample at T = 5.25 K. A 
Fourier spectrum analysis of these data shows that there 
are three maxima corresponding to Hj — 12.8 G, 9.9 G 
and 4.3 G which are relatively close to the matching con- 
ditions Hj/ Hi = 1.5, 1 and 1/2, respectively. Apparently 
these periodic jumps appear close to the boundary H* (T) 
with a jump size larger than that observed at lower tem- 
peratures and fields. 

These a priori unexpected commensurability effects at 
low temperatures represent a striking observation since 
just above the transition T,n p (H) there are no special 
features indicating the commensurability effects. The ab- 
sence of matching features above Tdi p (H) is an indication 
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FIG. 6: The main panel shows the decreasing branches of de- 
magnetization loops recorded at three different temperatures 
for the BH100 sample. In the flux-jumps regime the average 
magnetization tends to a common value. A full hysteresys 
loop is shown in the left inset. The right inset shows a zoom- 
in of the flux jump regime where quasi-periodic jumps appear. 
The units used in the vertical axes of the insets is the same 
as the one used in the main panel. 



that vortices are no longer regularly distributed. How- 
ever, the influence of the periodic pinning array might 
still play an important role as it favors an easier motion 
of vortices along the rows of the array than at any other 
orientation.- 1 This picture is consistent with the MO 
measurements of Ref.|l3| showing preferential penetra- 
tion along the principal axes of the pinning array. More- 
over, recent MO images^ performed on YBa2Cu307_a; 
thin films with a square array of holes at very low tem- 
peratures (T = 4.5 K) where no matching effects are 
observed, showed that the penetration of the flux front is 
highly anisotropic as a consequence of a strong guidance 
of vortices by the underlying pinning structure. 

Currently it is still a puzzling question why the pe- 
riod of the jumps coincides with multiples of H\. In 
Ref.^H the authors postulate that at low temperatures, 
regular vortex structures may appear in a flux depleted 
region near the border of the sample as a result of geomet- 
ric barriers (GB). As field is increased, this vortex-poor 
contour region is progressively filled with vortices which 
form a highly stable vortex pattern similar to those nor- 
mally observed close to T c . Above a certain threshold 
(i.e. Hj = nHi) the vortex distribution in the border be- 
comes unstable and is pushed towards the center of the 
sample by the screening currents thus triggering a flux 
avalanche. A first step in order to elucidate the origin of 
the observed quasi-periodic jumps is to discern whether 
GB are actually present in these samples. 

Some evidence of the presence of GB comes from the 



fact that matching features obtained at high tempera- 
tures are systematically better resolved when decreasing 
field than increasing field as a consequence of a delayed 
entrance of vortices in the sample but without affecting 
their exit. However, in most cases the observed differ- 
ence is almost imperceptible thus suggesting that the 
ubiquitous GB are not relevant. This result is consis- 
tent with the highly symmetric dc-magnetization loops 
observed for all temperatures studied as shown for ex- 
ample in the left upper inset of FigEI I n other words, if 
GB were indeed the responsible for the observed periodic 
jumps, this periodicity should be absent in the decreas- 
ing branch of the loop, in contrast to the general obser- 
vation. Additionally, in a recent work we have shown 
that the ac-susceptibility response can be accurately de- 
scribed within the simple Bean critical state model for 
the perpendicular geometry without invoking the exis- 
tence of GB. Indeed, in a system dominated by the GB, 
currents are constrained to flow along the sample's edge 
and therefore the sample's response could be modeled as 
that produced by a ring-shaped sample.— This particu- 
lar geometry will lead to a dome-like Cole-Cole plot more 
symmetric around x' — 0-5 than that shown for example 
in Fig. □ 

We can gain further insight in this particular issue by 
measuring the third harmonic susceptibility X3(T 7 H, h). 
The advantage of measuring this component lies in the 
asymmetry of the penetration-exit process which leads 
to a different shape of the minor loops traced out dur- 
ing an ac-cycle for the case of GB and bulk pinning. 33 
Fig. [7] shows one of these measurements for the AD65 
sample at H = 500 G and h = 0.3 G. Although we have 
also collected X3(T) data in the temperature range 5 K 
< T < 7.5 K, at several dc fields (0 G < H < 1 kG) and 
ac-drives (0.03 G < h < 3 G), for all these conditions 
the curves are very much alike. In the same figure we 
have included the temperature dependence of the third 
harmonic components X 3 (T) and X 3 '(T) calculated using 
the Bean model for a disc shaped sample with field ap- 
plied perpendicular to the plane of the diskiS2i2i As we 
can see, the theoretical curves reproduce qualitatively the 
main features of the measured third harmonic. In con- 
trast to that, these features cannot be accounted for by 
using the expression for a ring (similar to that obtained 
for edge barriers). This behavior becomes more evident 
in the inset of Fig. [7| where a Cole-Cole plot of the third 
harmonic using the same data shown in the main panel is 
presented together with the theoretical curve for the crit- 
ical state in a disk shaped sample within the Bean model 
scenario. We also included in this figure the curve for a 
ring rescaled by a factor xo to fit with the experimental 
data and accounting for the expected amplitude depen- 
dence in the case of GB (strictly this is not valid but it 
rescues the basic heart-shaped curve obtained in the case 
of GB which allows one to make a fast comparison). 

In the inset we clearly see that as for the disk geom- 
etry, the measured curve occupies the quadrants II, III 
and IV. In contrast to that, in the case of GB dominating 
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FIG. 7: Main panel: Temperature dependence of the real 
(filled circles) and imaginary (open circles) components of 
the third harmonic susceptibility for the AD65 sample at 
H = 500 G and h — 0.3 G. The theoretical expectations for a 
disk in a critical state described by the Bean model are shown. 
In the inset the theoretical curve (solid line) and experimen- 
tal data (filled circles) already shown in the main panel are 
represented in a Cole- Cole plot together with the curve for a 
ring shaped sample (dotted line). 



behavior, the first quadrant (corresponding to low tem- 
peratures) should also be covered. This feature remains 
for all amplitudes and dc fields studied, thus suggesting 
that the contribution of GB to the total response is very 
weak. On the other hand, the observed difference (mainly 
in x'-i{T)) between the Bean prediction and our experi- 
mental results, might be accounted for by assuming the 
more realistic case of field dependent critical current as 
in the Kim modeh 35 ! 36 

An alternative explanation for the observed periodic- 
ity in the flux jumps was recently proposed by Hebert 
et ali^. According to these authors, at low tempera- 
tures the flux profile corresponds to a multiterrace critical 
state composed by steps of constant B and zero critical 
current connected by abrupt changes in the flux density 
where the current is higher than that obtained from the 
average slope of the flux profiled In this scenario, the 
sudden penetration of a new terrace induces the move- 



ment of the internal terraces giving rise to jumps in the 
magnetization. At very low temperatures the concept of 
terraced state is lost as a more disordered flux distribu- 
tion appears. On the other hand since the avalanches are 
triggered by the local slope rather than the average slope 
of the flux profile, this model satisfactorily accounts for 
the observed increase of the H — T region where flux- 
jumps are detected in the patterned samples. However, 
a complete explanation of this puzzling observation re- 
mains elusive at the moment, and clearly further experi- 
mental and theoretical studies are needed. 

In the final stage of preparation of this manuscript we 
learnt about a recent report by Zhukov et ali 3 ^ showing 
several common results with the present work. 



IV. CONCLUSION 

We have shown that the low temperature magnetic be- 
havior of Pb films is dominated by vortex avalanches re- 
gardless the details of the pinning array. Although the 
size of the flux jumps is strongly reduced by introduc- 
ing a periodic pinning array, the region of the H — T 
plane where flux-jumps occurs is enlarged for patterned 
samples. This drawback can limit at low temeratures the 
applicability of these arrays to reduce the noise in SQUID 
systems as proposed recently; 10 At low temperatures and 
fields an almost T- and ii-independent magnetization 
is found. Although there is no clear clue to interpret 
this behavior, it might be related with the channeling of 
vortices by the predefined dendritic structures. Finally, 
we have demonstrated that the observed quasi-periodic 
jumps are unlikely to be originated in geometric barriers. 
Instead, a multiterrace critical state could satisfactorily 
account for this effect. 
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